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Continuing previous work, a model independent analysis of the solar neutrino anomaly is 
performed in terms of neutrino oscillations, allowing a comparison with the predictions 
of the Standard Solar Model. SMA and LMA solutions emerge also in this case, although 
somewhat different from the standard ones. The significance of the NC/CC double ratio 
measurable in SNO is illustrated in this context. 
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1 Introduction and motivations 

Flavour neutrino oscillations continue to be a pretty con- 
troversial matter. It is fair to recall that, so far, no direct 
signal of them, like neutrino appearance or explicit os- 
cillation patterns, have been observed cither in solar or 
in atmospheric neutrinos. The up/down asymmetry of 
the flux of the atmospheric and neutrinos gives, 
however, an indisputable evidence for the presence of 
an atmospheric neutrino anomaly. No equally clear ev- 
idence has been found, on the other hand, for the solar 
neutrino anomaly ^ ^, IJ] , since its standard interpre- 
tation relics on a combination of many different exper- 
imental and theoretical ingredients. Furthermore, the 
LNSD result still awaits for an independent confirma- 
tion. 

On the solar neutrinos, which are the subject of this 
paper, quite different attitudes can be taken, depending 
on the weight one gives to the input of the Standard 
Solar Model (SSM) §. On one side, it does not look 
reasonable to consider the solar neutrino anomaly as an 
artifact due to a large unknown error in solar models or 
in solar neutrino experiments. The other extreme atti- 
tude is to assume that all the ingredients of the analy- 
sis are correct, thus obtaining a rather precise determi- 



nation of the neutrino oscillation parameters. As well 
known, the best fits of the solar neutrino deficit in this 
framework are given by few peculiar energy-dependent 
survival probabilities. 

The truth is that unfortunately, so far, SupcrKamio- 
kandc (SK) has not found any evidence for a distortion 
of the energy spectrum, nor for Earth regeneration ef- 
fects, nor for seasonal variations of the neutrino flux. 
Furthermore, the most recent SK data |^ worsen the 
equality of the best fit, with the net result that the new 
best fit regions now include values of the oscillation pa- 
rameters previously discarded on the basis of the sole 
neutrino rates. Recent analyses found that all the dis- 
tinct best fit solutions have a high goodness-of-fit prob- 
ability j^. However, at least in part, this is just a re- 
flection of having fitted the few really problematic data 
together with many other 'degrees of freedom' that have 
not much to do with the problem. To really judge the 
quality of the fit one should perform a more complete 
statistical analysis or rewrite the data in terms of a min- 
imal set of 'optimal' observable^. 

In view of this situation, we find it useful to come 
back to an analysis which has minimal dependence upon 
the SSM inputs. This is the purpose of this paper, 



*The addendum at page ^ (seetion 5) is not present in the 
pubhshed version of this paper. 

tThe addendum at page ^| (seetion 6) is not present in the 
pubhshed version of this paper. 



*A similar comment can be done for atmospheric neutrinos. 
It is hard to judge if the — ♦ i^storilo interpretation gives an 
acceptable fit by looking only at the minimal of ^ global fit 
that includes electron data, low energy data and too many zenith 
angle bins. 
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continuing previous work along similar lines. From an 
experimental point of view, the main new information 
comes from the SK measurements, mentioned above, of 
the energy spectrum and of day/night or seasonal varia- 
tions of the neutrino flux. Their interpretation has little 
to do with the theoretical input of the SSM. 

The SSM independent analysis is performed in sec- 
tion 2. In section 3 we discuss its implications for new 
solar experiments. Conclusions are drawn in section 4. 
In appendix A we describe the details of the computa- 
tion. In appendix B we discuss how KamLand and neu- 
trino factories can test a high value of Amjg ^ 10"'^ eV^, 
allowed by solar data in presence of an undetected sys- 
tematic error in the Chlorine experiment. 

Fitting the solar, atmospheric and LSND anomalies 
with neutrino oscillations consistently with all boimds 
would require more than 3 neutrinos and peculiar mod- 
els. We limit ourselves to oscillations between the 3 SM 
active neutrinos and we await for a confirmation of the 
LSND result disregarded in the following. We use 
the same notations as in B. The three neutrino masses 



rrii are ordered such that Ar 



rameterized as 



'23 



> Ar 



> where 



The neutrino mixing matrix is pa- 



V - i?23(^^23) diag(l, e*"^, l)Rl3{dl3)Rl2{0l2) (1) 

where Rij (Oij ) represents a rotation by 9ij in the ij 
plane, < % < tt, and is a CP-violating phase. With 
these notations, 623 and Arn^^, ~ Am^g are relevant 
to the atmospheric neutrino anomaly, 612 and Amjj to 
the solar anomaly, while 6*13 can affect both solar and 
atmospheric neutrinos. 



2 SSM independent analysis 

One can perform a useful almost SSM independent anal- 
ysis H by just treating the overall and ^Be fluxes 
as unknown quantities, to be extracted from the data. 
Here we briefly recall how this procedure is justified 
(see 1^ for explanations and references). First, it is safe 
to use the standard spectral functions for the energy 
distributions of the single components 

i = {pp, pep, ^Be, ^^N, ^^O, ^^F, hep} 

of the solar neutrino flux, while the total flux of each 
component is regarded as unknown. Second, it is safe 
to set to their standard values the ratios ^l^n/^isq and 
$pep/$pp, to neglect ^''F neutrinos and to consider hep 
neutrinos only when computing the upper tail of the en- 
ergy spectrum of recoil electrons in SK. Although to a 
somewhat lesser extent, it is also safe to set $i3n/$7B(, 
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Figure 1: Best fit values of the neutrino solar fluxes, as 
obtained by this analysis, compared with SSM theoretical 
predictions. 



to its standard value. Finally the solar luminosity con- 
straint allows to express the pp flux in terms of the re- 
maining free parameters ^r^^ and <I>sb. 

This kind of analysis is useful because solar neu- 
trino rates have been measured by three different kinds 
of experiments. For any given oscillation pattern, each 
measured rate gives an allowed band in the ($83, ^^Be) 
plane (few examples arc shown in fig. Requiring 
a crossing of all the three experimental bands selects 
specific oscillation patterns. In this way one converts 
experimental data into informations on the oscillation 
parameters and on the neutrino fluxes $83 and ^tbc- 
This kind of analysis will become more powerful when 
the SNO and Borexino experiments will present their 
data. Already now the results are much more restric- 
tive than two years ago Q . SK and Gallium experiments 
have measured more precisely their fluxes and the new 
SK data now exclude in a SSM independent way a large 
part of the oscillation parameter space where MSW ef- 
fects are large. 

Furthermore, the CHOOZ bound on Pf, disappear- 
ance 1^ now holds for all values of the mass splitting 
Am^3 allowed by the SK atmospheric data. This was 
not the case one year ago, and implies that ^13 is small, 
6*13 < 15° at 95% C.L. Therefore 6*13 can only have a mi- 
nor impact on solar neutrino experiments. Unless oth- 
erwise indicated we will assume that ^13 = 0. 

The best- fit values of the neutrino fluxes $83 and 
$732 are shown in flg. |l|. The regions delimited by con- 
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Figure 2: Values of the solar neutrino fluxes (^"837 *J'~Bc) rneasured by the Chlorine experiment (continuous lines), 
the Gallium experiment (dashed lines) and by the SuperKamiokande experiment (long dashed lines) assuming four 
neutrino oscillation schemes: • the standard best-fit LMA point in fig. ^a; • the solar-model independent best-fit 
LMA point in fig. ^6; • the best-fit SMA point in fig. ^c; • the best-fit LOW point in fig. ^d. All four plots have 
the flux o/^B neutrinos in 10^ cm^^s""'^ on the horizontal axis and the flux o/^Bc neutrinos in 10^cm~^s~^ on the 
vertical axis. The ellipse is the Standard Solar Model prediction. All errors correspond to one standard deviation. 



tinuous (dashed) lines give the best fit at 90% (99%) 
C.L.[|The ellipses represent the 90% and 99% C.L. SSM 
prediction for these fluxes ||]: 

^bbIssm = 5.15 {Ito.ll) ■ cm\ (2a) 
$-Be|ssM = 4.8(1 ± 0.09) -lOVcm^s (2b) 

A standard analysis would include these theoretical con- 
straints in the x^, forcing <i>SB and $730 to be close to 
the SSM predictions. 

Fig. |l| shows that the best fit regions are neither far 
from the SSM predictions of eq. (||) nor peaked around 
them. This reflects the fact that oscillation patterns 
that gave the best standard fits of the measured neutrino 
rates are now disfavoured by the SK SSM-independent 
data. Basically there are two distinct best-fit regions in 

fig.0: 

• A region with (f>SB > 5 • 10^/ cm^s and $7Bo < 
5 • 10^/ cm^s produced by values of the mixing pa- 
rameters around the LMA solution. Fig. |^ shows 
how a perfect crossing of the three experimental 
bands occur around $8b ~ 7.5 • 10^/ cm^s and 
cf)7Bc ~ 3 • 10^/ cm^s. This crossing is obtained for 
Amfa = 4 • 10"^ eV^ and 6*12 = 0.42. The stan- 
dard analysis requires a crossing centered around 

tThc contour lines are drawn at Ax^ levels that correspond to 
90% and 99% C.L., if one converts values of Ax^ into "best fit 
probabilities" p using the standard expressions valid for a gaussian 
probability distribution. This is not a good approximation since, 
as frequently happens in solar neutrino fits, one finds few separate 
best-fit solutions, while a gaussian would have only one peak. A 
proper treatment would shift the values of 1 — p by relative 0(1) 
factors. A comparable shift would arise if we performed an exact 
marginalization of the joint probability distribution with respect 
to the oscillation parameters. We neglect such corrections, since 
they are comparable to the uncertainties of Bayesian inference 
arising from the need of choosing some prior distribution function. 



the SSM prediction: the best fit is obtained for 
a slightly larger values of 612 and gives the worse 
crossing shown in fig. |^. 

• A region with (f>SB £ [2. 5... 4] • 10^/cm^s pro- 
duced by values of the mixing parameters around 
the SMA solution. The best crossing, obtained for 
9i2 = 0.025 and Am^2 = 0.5 • lO'^ eV^ is shown 
in fig. ||c. It also gives the best standard fit. The 
previous best standard fit had larger 0i2 = 0.04 
and gave a crossing perfectly centered on the SSM 
prediction (see fig. Ic in [^), but is incompatible 
with the day/night and spectral SK data. 

Oscillation patterns around the LOW region (i.e. the 
one with large 612 and Am^ < 10^^ eV^) give a modest 
fit with $8B w (4 ^ 5) • 10^/ cm^s (see fig. |d). 

Before going on, it is useful to consider the region 
around <i>SB ~ 3 • 10^/ cm^s and <1>7bo ~ 0. This re- 
gion appears due to a unfortunate weakness of our SSM- 
independent analysis: assuming no oscillations, the three 
bands perfectly cross at $8b ~ 3 • 10^/ cm s and *1*7bc 
slightly negative. Therefore the no-oscillation case can- 
not be excluded at a high confidence level and vari- 
ous oscillations patterns not much different from the 
no-oscillation case provide acceptable fits. We consider 
such crossings as unfortunate accidents. Before fitting 
the mixing angles, we exclude by hand such cases by im- 
posing $8B > and $7Bc > 1-5 10^/cm^s rather than 
$7Be > 0. This does not conflict much with our purpose 
of performing a SSM-independent analysis, since very 
low values of the ''Be flux are unphysical, as the Boron 
neutrinos, seen in SK, originate from the Berillium ones 
to a large extent. 

Fig. shows the flt in the usual plane of the mixing 
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Figure 3: Best fit value of the oscillation parameters (a) standard fit (h) solar-model-independent fit. 



parameters 9i2 and Amfj: the standard fit is sliown 
in fig. I^a and the the solar-model independent fit in 
fig. ^D. The best fit regions are not entirely restricted 
to 912 < 7r/4 jn). 

Concerning the standard fit, we mention an impor- 
tant detail not immediately apparent from the figure. 
Like other standard fits our fit still contains a SMA 
region, even if the SK spectrum and day/night data have 
excluded the 'old SMA' region with larger 9i2. The 'old 
SMA' gave such a good standard fit of solar rates that 
values of 9i2 previously regarded as 'too small' now give 
an acceptable standard fit of solar rates. This is why we 
obtain a new best-fit SMA region, more SMA than the 
old one (see fig. ^). Fig. ^ explains why such smaller 
values of 6*12 were discarded in old standard fits but not 
in old SSM-independent fits: they give a good crossing 
of the three experimental bands, but at a value of the 
Boron flux smaller than the one predicted by the SSM. 
This means that the SK spectrum and day/night data 
are not a problem for the solar-model independent SMA 
region]^. Non-zero values of 9i^ just below the CHOOZ 
bound slightly shift the crossing point towards higher 
Boron fluxes, and therefore slightly improve the quality 
of the standard fit in the SMA region. 

Presently the LOW solution gives a better standard 
fit than the SMA solution Q. Fig. ||d shows how the 
three experimental bands cross in the case of the 'best 

-fThis discussion implies that the sterile neutrino interpretation 
of the solar neutrino anomaly is disfavoured by SK only if one 
imposes the SSM value of the Boron flux. 



standard fit' LOW solution, 9i2 = 0.66 and Arnfj = 
0.8 10^^ eV^. The crossing is not good, but roughly 
centered on the SSM prediction. A solar-model inde- 
pendent analysis does not reward this property. The 
best SSM-independent fit in the LOW region has larger 
9i2 and lower ^tbc than in the standard fit. 

The band corresponding to the Ga experiment in 
fig. ^ (the almost horizontal one) is not unacceptably 
high because Earth-regeneration effects strongly affect 
neutrinos with energies w MeV(Am74 IQ-'^cV^). 
Unfortunately radiochemical experiments, which detect 
such neutrinos, cannot study day/night effects. Earth- 
regeneration gives a < 10% seasonal variation of the 
capture rate in GNO, since at Gran Sasso nights are 
longer in winter than in summer |12 . Gallex does not 



see such an effect. Present data from all Gallium exper- 
iments could be sensitive to a 10% seasonal variation. 



3 Expectations for SNO 

The fact that the SMA solution has migrated toward 
smaller values of 9i2 , previously considered only in SSM- 
independent analyses |^ , has significant implications for 
the SNO experiment. Previous studies of the signifi- 
cance of the SNO experiment found that only a 
global fit of various SNO precision observables could 
eventually discriminate between the SMA and the LMA 
solutions. 

On the contrary, we think that it is quite possible at 
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Figure 4: f 117 how the SK hounds on spectral and day /night effects (continuous line) shift the best- fit SMA 
region: the dashed area shows the old SMA, the gray area the new SMA. All lines correspond to best fit-bounds at 
90% C.L. Fig.s ^&,c show the probability distributions (see text) for the double NC/CC ratio in the new (fig. 
and old (fig. Qcj SMA and LMA best fit regions. 



SNO to discriminate between the new LMA and the new 
SMA solution| or even find evidence for oscillation pat- 
terns suggested by our solar-model independent analysis 
by making use of the NC/CC double ratio 

(NC rate)/(CC rate) 

r = . 

(NC rate)/(CC rate)no oscillation 

r is an interesting observable because dominant theo- 
retical uncertainties cancel out when taking the double 
ratio, and because expected oscillation effects can be 
very large so that one does not need a very precise de- 
termination. 

In order to perform a quantitative analysis, we com- 
pute the values of r and of the x^-probability 

p = exp(-AxV2) 

for a grid of oscillation patterns in the SMA and in the 
LMA regions. The Ax^ is computed with respect to the 
local LMA or SMA minimal so that p = 1 in the 
best- fit SMA point and in the best-fit LMA point. We 
assume an energy threshold > 5MeV on the recoil 
electrons originating from CC interactions Ved —> ppe, 
but our final results do not depend on this choice. 

In fig. ^,c we plot p{r), the maximal value of p at 
which any value of r can be reached in the LMA and 
in the SMA region. The x^ is computed performing a 
standard analysis. In fig. ^ we have included the most 
recent SK data, while in fig. ^ we have not included 
them. Fig. ^ shows that there is now a neat separation 
between SMA and LMA predictions: measuring r < 2 
or r > 2 would have clear implications. A measurement 
of r can also provide a signal for non standard solutions. 

§This observation was also made in a recent paper |l4|. 



Dividing the possible values of r in 5 distinct ranges, we 
can summarize the situation in the following way: 

1. values of r very close to 1 are allowed in the non 
standard part of the SSM-independent SMA re- 
gion (with smaller 6*12 and low flux). Oscilla- 
tions into sterile neutrinos would also give r ~ 1. 

2. 1 < ?■ < 2 is allowed in the standard or SSM- 
independent SMA region. 

3. r very close to 2 is the value predicted by the 
non standard solution with high ATO22 ^^'^ ^12 ~ 
7r/4 (8[ |l^, allowed in presence of an undetected 
systematic error in the Chlorine experiment. 

4. 2 < r < 4 is allowed in the standard or SSM- 
independent LMA region; 

5. 4 < r < 5 is still allowed in the non standard part 
of the LMA region; 

For completeness it should be said that values of r be- 
tween 1.5 and 3 are expected also in the standard LOW 
region. 

The SNO experiment will improve also the experi- 
mental knowledge of the solar neutrino fluxes. The NC 
rate is not affected by oscillations between active neutri- 
nos, and therefore provides a measurement of the Boron 
flux. It is expected to have a < 10% systematic error, 
mainly due to the uncertainty in the detection cross sec- 
tion jl^ . Due to the large spread between the values of 
the Boron flux required by the different oscillation pat- 
terns (see fig. |l|) this measurement will also have a sig- 
nificant impact. In particular the SMA solution requires 
low values of the Boron flux. 
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The Borexino experiment will be mainly sensitive 
to the Berillium component of the solar neutrino flux. 
Therefore its data will be represented by one quasi- 
horizontal band in fig.s|^, at a level dependent on the ac- 
tual oscillation pattern. Although the main features can 
be seen already from fig.s and a true understanding 
of the impact of Borexino data on a SSM-independent 
analysis will require a combined fit of the oscillation pa- 
rameters and of the neutrino fluxes. 

The KamLand reactor experiment will measure 
accurately the oscillation parameters, if they lie in the 
LMA region (see fig. ^). In this case, the solar neutrino 
data will give the ^B and ^Be solar fluxes. In partic- 
ular, the Borexino data will be crucial for an accurate 
determination of the ^Bc flux. 

4 Conclusions 

The SK measurements of the energy spectrum and of 
day/night and seasonal variations of the neutrino flux 
have not realized, so far, any "smoking gun" in the study 
of the solar neutrino problem. Nevertheless these mea- 
surements provide significant information, since inde- 
pendent from theoretical models. Their use, combined 
with a proper treatment of all the different rate mea- 
surements allows an almost direct determination of the 
preferred values of the ®B and ^Be solar neutrino fluxes. 
In turn these values can be compared with the SSM ex- 
pectations. 

This comparison at present is encouraging but far 
from conclusive, as illustrated in fig |l|. In particular, 
also in view of fig. ^ it makes it clear how premature it 
is to select one specific pattern of neutrino oscillations to 
explain the solar neutrinos. Nevertheless, data expected 
in the near future especially from SNO or from Borex- 
ino can turn this comparison into a convincing proof of 
solar neutrino oscillations and can also provide, at the 
same time, an independent validation of the SSM from 
neutrino physics. 

Acknowledgments We thank Andrea Romanino for 
useful discussions. 

A Details of the computation 

The energy spectra for the independent components of 
the solar neutrino flux have been obtained from 
The neutrino production has been averaged for each 
flux component over the position in the sun as predicted 
in 1^, |l^ . This averaging does not give significant correc- 
tions. MSW oscillations inside the sun have been taken 



into account in the following way. The 3x3 density ma- 
trix PS for neutrinos exiting from the sun is computed 
using the Landau-Zener approximation with the level- 
crossing probability appropriate for an exponential den- 
sity profile [0 The density profile has been taken 
from and is quasi-exponential: small corrections to 
ps have been approximately included. Oscillation effects 
outside the sun are described by the evolution matrix U, 
so that at the detection point pE = UpsU^ ■ In particu- 
lar, earth regeneration effects have been computed nu- 
merically using the mantle-core approximation for the 
earth density profile. We have used the tree level Stan- 
dard Model expression for the neutrino/electron cross 
section at SK. The CC and NC cross sections at SNO 
have been taken from . The experimental resolution 
at SK and SNO has been included as suggested in . 

The total neutrino rates measured with the three 
kind of experimental techniques arc ^ ^ Q 

*ci|exp = (2.56 ± 0.22) 10-36s-i (3a) 
$Ga|exp = (74.7±5) lO-^V^ (3b) 
^skIcxp (2.40 ± 0.08) • 10^ cm-^s-i (3c) 

having combined systematic errors in quadrature with 
statistical errors. The SuperKamiokande experimental- 
ists give directly the value of the fiux they measure. 
The other experiments involve more uncertain neutrino 
cross sections and prefer to give the frequency of events 
measured per target atom in their detector. 

The solar model independent SK data included in 
the fit are: the energy spectrum of the recoil electrons 
(divided in 18 energy bins between 5.5 MeV and 15 MeV) 
and the total flux measured at SK during the day and 
during five night bins (defined according to the value of 
the cosine of the zenith angle) H, |^. The SK collab- 
oration can include in their fit data about the zenith 
angle variation of the recoil electron spectrum and ex- 
clude the old SMA at 95% C.L. With these unpubhshed 
data the standard and the SSM-independent SMA so- 
lution would be less attractive and fig. |l| would show a 
more neat separation into two distrinct regions. 

B Large Amf2 and nu-factories 

The standard interpretation of the solar neutrino anomaly 
is based on many experimental and theoretical ingredi- 
ents. We have discussed how the SSM predictions can be 
tested. There is one other ingredient that could be not 
solidly founded and that has a significant impact on the 
final result ^ |l^ . Only a single experiment, the Homes- 
take one, has detected neutrinos with the Chlorine tech- 
nique (with the other techniques, data come from two 
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water Cerenkov and two Gallium experiments). Fur- 
thermore Homestake is the only experiment that ob- 
serves a rate different than one half of the SSM pre- 
diction in absence of oscillations, therefore excluding an 
energy-independent survival probability Pee{E^) « 1/2. 
This important conclusion could be the result of an 
under-estimation of the systematic error, that according 
to the Chlorine collaboration is equal to the statis- 
tical error. It would be interesting to perform a direct 
calibration of the Chlorine detector ||2^ . 

Pee{Ey) ~ 1/2 can be obtained with 612 ~ 7r/4 and 
^"^12 ^ 10^^ eV^. This oscillation pattern has no prob- 
lems with the recent SuperKamiokande data so that, 
even in a standard analysis, it is no longer significantly 
worse that the new best fits. 



0.8 



0.6 



0.4 



0.2 



7 8 i^gS 2 3 4 5 6 7 8 ^q4 

Pathlength L in km 



The KamLand experiment will test the LMA 
region looking at disappearance of Ve reactor neutrinos. 
If 

A77i^2 fmd 0\2 he inside the LMA region, KamLand 
can accurately measure them 22 . If instead Am^2 ^ 2 • 
10"'' eV^, i>e oscillations are averaged so that a measure- 
ment of Ato^2 needs a good energy resolution, a precise 
knowdlege of the un-oscillated spectra, high statistics 
and low background. Assuming that these conditions 
can be satisfied, fig. ^ shows the accuracy at which Kam- 
Land can measure few values of Amfj and O12 (repre- 
sented by the dots) after three years of running (i.e. with 
2400 events if no oscillation occur) . If Am^2 is too large 
statistical fiuctuations often lead to discrete ambigui- 
ties in its determination. A reactor experiment with a 
shorter baseline would not have this problem. 

Here, we study the impact of a large Am^2 at a neu- 



trino factory 



Due to the high beam purity, a neu- 



trino factory will allow to study v^. — > and v,, 
oscillations down to small values of the oscillation prob- 
ability. Extensive studies |2^ have determined the 
optimal energy and pathlength that give the maximal 
sensitivity to a small ^13. If ^13 = 'solar' oscillations 
give effects oc (Ato22)^ that can be seen at a neutrino 
factory if Ato?2 ^2.5 • 10"" eV^ With a non-zero 6*13 
and even for values of ^m\2 inside the LMA region, 
'solar' effects oc Am^2 have a significant impact on 6'i3 
measurements at a neutrino factory p^ . 

The most promising observable for discovering 
oscillations is given by fj." appearance at a relatively 
short baseline L w 700km. The number N{ii^) of fi^ 
events produced by both 'solar' and 'atmospheric' oscil- 
lations can be approximated by treating 6*13 and A77i^2 



effects as perturbations. This give^ P{i^e 



Figure 5: Suppression 0/ Am^2 oscillations due to mat- 
ter effects as function of the pathlength L 
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(4) 



where iV^ is the number of decays occurring in the 
straight section of the storage ring pointing towards the 
detector, ~ (20 50) GeV is the energy, A^kt is 
the size of the detector in kilo tons, e is the efficiency for 
the detection of /x^ . Using the parameterization (0) 
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where A = 2V2NeGFE^, A' 



A, a 



cos 6j 



and Sij = sin dij. If the phase factors are large one 
should average N(fj,~) over the neutrino spectrum, oth- 
erwise one can set E^ « E^. For small L (in practice for 
L < 700 km) reduces to the e/x element of the neu- 
trino squared mass matrix. Assuming a short baseline, 
Am^2 < Ami3, 6*12 « 6*23 « 7r/4 and 6*13 < 1 
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we clearly see in an analytical way how significant Am^2 
oscillations can be. An analogous approximation holds 



for Vf, 



signals. 



At baselines L > 10'^ km matter effects become signif- 



t Using the formula e^^+' 



leant. Using eq. (||), the number of events due to At 
'i^+J^ e"^^' -e-e^^'dx+Oie"^)). oscillations Only can be written as its value at L = 
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Figure 6: Few fits at 90% and 99% C.L. of simulated 
KamLand data after three years of running. 



multiplied by the function 

sin^ X 
V{L) = < 1, 



which exact numerical value is plotted in fig. ^ There- 
fore a short baseline L w 700 km is optimal for discov- 
ering Am^ 



1^2 effects. 



By performing a global fit of simulated nu-factory 
data wc find that it will be difficult to distinguish Ato^2 
effects from 6'i3 effects at a good C.L. by comparing data 
taken at different pathlengths and/or different neutrino 
energies, as suggested in p3]. For certain values of the 



CP violating phase a comparison between 



and 



i^/j rates allows a better discrimination. Such 'pre- 
cision studies' are statistically significant only with a 
sufficient number of events. For example, observing few 
events only would not allow to say if they are due to a 
013 ~ 0.007 around its nominal sensitivity, or due to a 
Amfa « 3 lO"'* eV^. In conclusion, if Am^ > 2 10-4 eV^ 
so that KamLand cannot measure it, an accurate mea- 
surement of Am^ cannot even be obtained with a neu- 
trino factory: a new reactor experiment with intermedi- 
ate baseline ~ 10 km would be necessary. 

If KamLand will give a precise measurement of Am^j 
free from discrete ambiguities (this could not be the case 
if Am22 > 2 • 10-4 eV^), by comb ining KamLand data 
with nu-factory data one can usually obtain a satisfac- 
tory fit of 013 and of the CP violating phase. 
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Figure 7: Fig. || updated including the SNO CC re- 
sult: best fit, at 90% CL (continuous line) and 99% CL 
(dashed line), of the neutrino solar fluxes compared with 
SSM theoretical predictions. 

5 Addendum: SNO CC results 

In this addendum, we update our results by adding the 
first SNO data. The SK cohaboration published the 
day + night spectral data Q as 19 + 19 energy bins, 
so that we can now include this information in the fit. 
We also exphcitly included the CHOOZ data (l0|, that 
disfavour values of Ato^j above 0.7 10"'^ eV^ for large 
mixing 9i2 ~ 7r/4. Finally, we improved the numerical 
accuracy of our computation, and extended it to include 
vacuum oscillations 

As pointed out in section 3, the measurement of the 
NC/CC ratio alone was expected to discriminate be- 
tween LMA and SMA. The measured value happens to 
fall into case 4. (cases 3. and 5. are not significantly dis- 
favoured). Therefore the SMA region is now strongly 
disfavoured, even from a solar-model-independent point 
of view. Furthermore, the solar-model- independent LMA 
region (case 5.) no longer gives a best significantly 
lower than the standard fit. These results are confirmed 
by an updated analysis, as shown in fig. |^ and ^. 

In fig. 1^ we show the updated fit of the solar neutrino 
fluxes. Values of the flux detectably different from 
the SSM prediction no longer give good flts. On the 
contrary, a discrepancy between the ''Be flux and its 
SSM prediction could still have significant effects. The 
best-fit region at 90% CL in fig. ^ is composed by two 
disjoint regions. The one with smaller ffuxes is obtained 



from vacuum oscillations (but can also be obtained with 
LMA, LOW, SMA oscillations with a worse CL). The 
one with larger fluxes is obtained from LMA oscillations. 

In flg. H we show the updated global fits for the os- 
cillation parameters: flg. |8^ shows a standard fit, while 
flg. |8|d is the solar model independent flt described in 
section 2. In flg. ||c we perform a standard fit, but drop- 
ping the uncalibrated Chlorine rate from the fit. The 
results of fig.s ^,b,c are quite similar: more or less ac- 
ceptable fits can be obtained for a large range of Am^ 
and large mixing angle, while the SMA solution is al- 
ways disfavoured. 

Finally, we mention another important aspect of so- 
lar model independent analyses. It is sometimes said 
that it is useless to measure the Gallium rate with an 
error much smaller than the solar model uncertainty. 
As explained in section 2 and illustrated in flg.s ||, this 
is not true. Fig.s ^ in fact shows that the amount of 
information that can be extracted in a solar model in- 
dependent way from the solar rates is today limited by 
the accuracy of the Chlorine experiment. This limita- 
tion will disappear when Borexino will measure the ^Bc 
flux. Fig.s H show that solar model independent con- 
siderations already now give useful informations on the 
oscillation parameters. 
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Figure 8: Global fits of solar data (updated including SNO CC data) at 90, 99% CL in the (tan^ 6*12, Am^2/ eV^) 
plane: (a) standard fit; (b) solar model independent fit; (c) standard fit dropping the uncalibrated Chlorine data. 
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Figure 9: Global fits of solar data (updated including SNO NC data) at 90, 99, 99.9% CL in the (tan^ ^12, Am^2/cV^) 
plane: (a) standard fit; (b) solar model independent fit; (c) standard fit dropping the uncalibrated Ghlorine data. 



6 Addendum: SNO NC results 

In this addendum, we update our results by adding the 
latest SAGE, GNO and SK data and the SNO day/night 
energy spectrum As explained in the SNO paper, the 
SNO spectrum allows to extract its NC, CC and ES con- 
tributions. For example, assuming energy-independent 
oscillations between active neutrinos, our reanalysis of 
SNO data gives 

$Nc = (5.2±0.5) lOVcm^s, (5a) 

$cc= fee$8B = (1.76 ± 0.08) 10^ cm^s. (5b) 

(the errors are somewhat anti-correlated) .[] 

^Our reanalysis slightly differs from the corresponding SNO 
analysis in jlj. SNO extracts from all data (energy and zenith- 
angle spectra) the CC, NC and ES components. We instead ex- 
tract the CC, NC components from the energy spectrum, assum- 



As expected, new data have a significant impact on 
the oscillation fit Q. In fig. ^ we show the updated 
global fits for the oscillation parameters: fig. ^ shows a 
standard fit, while fig. ^ is the solar model independent 
fit described in section 2. In fig. ^ we perform a stan- 
dard fit, but dropping the uncalibrated Chlorine rate 
from the fit. The results of fig.s |^a;b,c are quite similar: 
more or less acceptable fits can be obtained for a large 
range of Am^ and for large mixing angle, while the SMA 
solution is always incompatible with data at about 5a 
level. In all analyses LMA oscillations give a better fit 
than LOW and (Q)VO oscillations. In fig.s ^,b,c we 

ing the standard relation <1?es ~ ^CC + O.IS^nc (the ES rate 
have also been accurately measured by SK). We take into account 
systematic errors and backgrounds as described by SNO [0, as 
well as the theoretical uncertainty on the Boron energy spectrum 
(that also affects other solar neutrino experiments). See for a 
useful discussion of these issues. 
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Figure 10: Fig.s |^ updated including the SNO NC 
result: best fit, at 90% CL (continuous line) and 99% 
CL ( dashed line ), of the neutrino solar fluxes compared 
with SSM theoretical predictions. 

have not included the CHOOZ bound, in order to show 
that Am^ > 10"'^ eV^ is now disfavoured by solar data. 

In fig. |l^ wc show the Boron and Beryllium fluxes, as 
extracted from the solar neutrino data without assuming 
solar model predictions for these fluxes, and assuming 
that the solar anomaly is due to i^e — > i^p,,T oscillations. 

Within this framework, SNO has measured the Boron 
flux through NC scattering, finding a value in agreement 
with solar model predictions. A much larger Boron flux 
was already excluded in our previous analysis (see fig- 
ures |l] and 1^) because it would need a small and energy- 
independent survival probability Pee(-E^i/)i that cannot 
be obtained by oscillations. 

In view of this experimental progress, one can study 
if new significant information, e.g. about the CNO and 
pp fluxes, can now be extracted from a more general 
analysis. Since the results are not very interesting, we 
just briefly describe their main features. Only two kind 
of experiments. Gallium and Chlorine, have measured 
low-energy neutrino fluxes Therefore, one can think of 
extracting the values of two of these fluxes at most. The 
Chlorine experiment is not very sensitive to low energy 
neutrinos: after subtracting the ^ 80% Boron contribu- 
tion, as measured by SNO via CC, the residual Chlorine 
rate is just about 2cr above zero. 

A first interesting question is: do solar neutrino data 
discriminate between CNO and ^Be neutrinos? We find 



that equally good fits are obtained in the extreme limits 
of vanishing CNO or ^Be flux. Therefore fig. |l^, where 
the Be flux is studied, provides all the significant infor- 
mation on Be and CNO fluxes, and implies that their 
sum is at most 2 times larger than solar model predic- 
tions. This means that the solar luminosity constraint 
fixes the pp flux to be around its solar model value. 

A second issue is: what is it possible to say on the pp 
flux, without imposing the luminosity constraint (and of 
course without assuming the solar model predictions for 
the other fluxes). The relatively more interesting result 
is that the pp flux cannot be zero, as can be seen by 
comparing the Gallium, Chlorine and SNO CC rates. 

Ref. 1^ finds that the pulls between solar model pre- 
dictions and their best-fit values (as obtained from a 
global oscillation fit that includes these predictions) are 
small. Note, however, that solar models are not con- 
firmed in this way. While a large pull would signal a 
problem, a small pull is obtained in two cases: when 
the experimental determination agrees with the predic- 
tion (this happens in the case of the Boron fiux) but 
also when the experimental error is much larger than 
the theoretical error (this roughly happens for all other 
fluxes). Looking at pulls only, one cannot discriminate 
these two extreme cases. In order to test predictions one 
must extract the predicted quantities (e.g. the B and Be 
fluxes in fig. |l^, more generically any other 'systemat- 
ics' relevant for solar data) by fitting the data without 
assuming the predictions under examination. 
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